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Abstract 

Given positive integers a,b,c and d such that c and d are coprime we show 
that the primes p = c(mod d) dividing a k + b k for some k > 1 have a 
natural density and explicitly compute this density. We demonstrate our 
results by considering some claims of Fermat that he made in a 1641 letter 
to Mersenne. 

Mathematics Subject Classification (2001). 11N37, 11R45. 



1 Introduction 

If S is a sequence of integers, then we say that an integer m divides the sequence 
if it divides at least one term of the sequence. The sequence {a k + b k }^ =l we will 
denote by S a> b. Several authors studied the problem of characterising (prime) 
divisors of the sequence S a ,b- Hasse [5] seems to have been the first to consider 
the Dirichlet density of prime divisors of such sequences. Later authors, e.g., 
Odoni [13] and Wiertelak strengthened the analytic aspects of his work, with 
the strongest result being due to Wiertelak [18J. In particular, Theorem 2 of 
Wiertelak [15] , in the formulation of [llj, yields the following corollary (recall 
that Li(x) = dt/\ogt denotes the logarithmic integral): 

Theorem 1 Let a and b be positive integers with a ^ b. Let N a ^(x) count the 
number of primes p < x that divide S a ^. Put r = a/b. Let A be the largest integer 
such that r = u 2> ~ , with u a rational number. Let L = Q(i/u). We have 

at ( \ xr \r-f \ , n ( ^(ioglogz) 4 

N a ,b{x) = d(r)Li(x) + O — - — = 

V log X 

where the implied constant may depend on a and b, and 8{r) is a positive rational 
number that is given in Table 0. 

Table 0: The value of 5(r) 



L 


A 


5(r) 


L ± Q(V2) 


A > 


2 x - A /3 


L = Q(V2) 


A = 


17/24 


L = Q(V2) 


A = 1 


5/12 


L = Q(V2) 


A > 2 


2~ A /3 
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Theorem [T] implies that if a and b are positive integers such that a ^ b, then 
asymptotically N a ^(x) ~ 5{r)x/ \ogx with 5{r) > 0. In particular, the set of 
prime divisors of the sequence {a k + b k }^L 1 has a positive natural density. 

In this paper we will establish, inspired by a letter from Fermat (see next 
section), a related result. 

Theorem 2 Let a, b, c, d be positive integers with (c, d) = 1 and assume that 
a b. Let r and A be as in the previous theorem. Let 

N a b(c,d)(x) := #{p < x : p\S ab , p = c(mod d)}. 



Then, for 



ah < log 2 / 3 x and d < — ^- , 

log log x 

N atb (c,d)(x) = 5 a:b (c, d)Li(x) +0 \ x ^°&^>& x 

V log ' x 

where S a ^(c, d) is a rational number that is given in Tables 1 to 6 and the implied 
constant is absolute. 

We have < 5 a ,&(c, d) < l/(p(d) by the prime number theorem for arithmetic 
progressions. In case S a ,b(c,d) = there could potentially be infinitely many 
primes p = c(mod d) dividing S a ^- However, using elementary arguments not 
going beyond quadratic reciprocity, one can show that there are at most finitely 
many primes p dividing S a> b in this case. Likewise if 5 a fi = l/(f(d), using ele- 
mentary arguments not going beyond quadratic reciprocity, one can show that in 
each case there are at most finitely many primes p = c(mod d) not dividing S^. 
For a more precise statement we refer to Theorem HI 

Inspection of the tables shows that we can always write ip(d)S a; b(c, d) = 7^3, 
for some non-negative integers c and m. 

Notations: 

As the tables for the density depend on some auxiliary parameters computed from 
a, b, c, d, some notations are needed to read them. We introduce these notations 
here and they will be maintained throughout this article. Given a, b and the 
modulus d, there is a unique table among the 6 from which one reads off the 
density. Put r = a/b = Tq, where ro is not a proper power of a rational number. 
Write h = 2 x h', d = 2 5 d', with h', d! odd. Put Vi{c— 1) = 7, where it is understood 
that 7 is larger than any number when c = 1. We denote the discriminant of the 
quadratic field Q(y/i) by D{t) and we put D(r ) = 2 5 °D'. We also write r = u/v 
and t = — r or Yli=i('^)Pi according as to whether uv is odd or uv = 2 ni=iP«- 
By d°° we denote the supernatural (Steinitz) number ripidP 00 - For each positive 
integer j > 1, we put N 3 = r 1 ^'" 1 , Q) and iVj = Q{( 2j , r 1 ' 23 , Q), where 

for any /, denotes any fixed primitive l-th root of unity. Finally, for j > 1, the 
intersection fields Kj := Q(( 2J ■, r 1 / 2 ^ 1 ) n Q(( d ) and K'- := Q(Q, r 1 ^) n 
will occur throughout our discussion. 



we have 
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Table 1 : Q(v^o) ^ Q(v^), D' \ d! 



X 


S 


0(d)M c > rf ) 


< s 


<7 


2 A+1-A 
3 


* 
* 


> 0, < min(A,7) 



3 

2 1-A 
3 


>7 


> 7 





<7 


> 7 


1 - 2 a -t 



Table 2 : Q(^) ^ Q( v / 2), Z>'K 5 < 5 



X 


5 




<j)(d)8 a;b (c, d) 


>5-l 


>0,< 7 


1 


3 






-1 


2 &-i-x 


* 





1 


3 






-1 




<5-l 


<7 


1 


1 3 






-1 


1 


> 5 


> 7 


* 





< 7 - 1 


>7 


1 


1 _ 2 a + 1 -t 






-1 


1 


>7 


> A 


* 






Table 3 : Q(^) ^ Q(v / 2), £>K and 8 > 5 



X 


5 




<t>{d)5 afi (c, d) 


<5-l 


< 7 


1 


2*+i-» 2 x + 2 + 6 - 26 
3 3 


<S-1 


<7 


-1 


2\+l-S 2 x +' 2 + / >-' M o 
3 3 


= 5-1 


<7 


1 


2 , 2 2d+L -' 2d o 

3 ' 3 


= 6-1 


<7 


-1 


2 2 2t >+ 1 -' 2d o 

3 h, — 


< 7 - 1 


> 7 


* 


1 - 2 a -t 


>7 


> A 


* 





> 5 


>7 


* 





<S -2 


> 0, < min(7, A) 


1 


2 A + 2 + d_2A 
3 3 


<5 -2 


> 0, < min(7, A) 


-1 


2<>-a 2 A + 2 + d - 2 *o 
3 3 


>5 -l 


>0,< 7 


1 




>5 -l 


> 0,< 7 


-1 




<S -2 





1 


2 1 — A 2 A + 3 — !2 d q 

3 3 


<5 -2 





-1 


2l — A <2\-\-'6 — 2d() 

3 2— t- 


>5 -l 





1 


2~ A 
3 


>S -1 





-1 


2-a 



3 



Table 4 : Q(^) = Q( v / 2), 5 < 2 



\ 

/\ 


u 


7 




U 


< 1 


^ o 


1 I J 24 


U 


2 


> 


11 / 1 o 

11/12 


U 


z 


1 


1/2 


1 


2 


1 





1 


< 1 


> 5 


5/12 


1 


2 


> s 


5/6 


> 2 


< 1 


> 6 


2" A /3 


> 2 


2 


> 5 


2^/3 


> 2 


2 


1 






Table 5 : Q(^o) = Q(V2), 5 > 3, A > 



A 


5 


7 


<p(d)5 atb (c, d) 


> 2 


3 


< 5 





> 5-1 


> 3 


> $ 




> 2,< <5- 1 


> 4 


> 5 


1 3 


>2,<7-2 


> 4 


< 5 


1 _ 2 A +!-7 


> max(2, 7 — 1) 


> 4 


< S 





1 


> 3 


> s 


-, 2'^ 

3 


1 


> 3 


1 





1 


> 3 


2 


1 


1 


> 3 


> 3,< 5 


1 - 2 2 ^ 



Table 6 : Q(y^) = Q(y/2), 5 > 3, A = 



7 


c(mod 8) 


4>(d)5 a ,b(c, d) 


> 5 


1 


1 2 2 -" 
3 


< 2 


±1 





< 2 


±3 


1 


> 3,< 5 


1 


1 - 2 x -t 



In the next section we reconsider a letter from Fermat and papers by 3 authors 
[H [2j [T7] in the light of Theorem [2j In Section 3 we prove Theorem [2J except for 
the fact that an expression for 5 a ^(c, d) in terms of data from algebraic number 
theory appears. In Sections 4-7 we evaluate this expression for 5 a i ,(c, <f). The 
outcome is recorded in Tables 1-6. This then completes the proof of Theorem 
|2j In Section 8 we determine the cases in which S a ,b(c, d) = 0, respectively 
8 a ,b{c, d) = l/f(d). In the final section we produce the results of some numerical 
experiments and show that they match well with what can be read from our 
tables. 
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2 On a letter of Fermat to Mersenne 



Fermat [U p. 220], cf. Dickson [3l p. 267], in a letter to Mersenne dated 15 June 
1641 stated that (p will always be used to denote primes): 

Conjecture 1 (Fermat, 1641) 

1) Ifp\S 3jl , thenp^ -l(mod 12). 

2) Ifp\S 3)1 , then p ^ +l(mod 12). 

3) If p\S^i, then p ^ — l(mod 10). 

4) Ifp\S 5)1 , then p ^ +l(mod 10). 

Pur r = a/b. For p \ ab there exists a smallest positive integer k such that 
r k = l(modp); this is ord p (r), the multiplicative order of r(modp). It is not 
difficult to see that if p \ ab, then p\S a ^ if and only if ord p (r) is even. If p\ab 
and p \ (a, b), then clearly p \ S a ^. (With (a, b) and [a, b] we denote the greatest 
common divisor, respectively lowest common multiple of a and b.) Using this 
observation and the law of quadratic reciprocity it is easy to see that the following 
holds: 

Proposition 1 Conjecture 1.1 of Fermat holds true. 

Proof. For p > 3 by the law of quadratic reciprocity we have (|)(|) = (— l)^ - - 
Suppose that p = — l(mod 12). It then follows that (5) = 1. By Euler's identity 

we then have 3^ = (-) = l(mod p). Since (p — l)/2 is the largest odd divisor 
of p — 1 it follows that ord p (3) is odd. This implies that p \ S 3> i. □ 

However, a computeralgebra computation learns that the remaining conjectures 
are all false. Counterexamples (in ascending order) are listed below : 
Counterexamples to: 



Conjecture 1.2: 37, 61, 73, 97, 157, 193, 241, 337, 349, 373, 397, 409, 457, ■ • • 
Conjecture 1.3: 41, 61, 241, 281, 421, 521, 601, 641, 661, 701, 761, 821, 881,- • 
Conjecture 1.4: 29, 89, 229, 349, 449, 509, 709, 769, 809, 929, 1009, 1049, • • • 



Sierpihski suggested that Conjecture 1.2 is false for infinitely many primes. This 
was proved by Schinzel [T6], who in the same paper showed that also Conjecture 
1.3 and Conjecture 1.4 are false for infinitely many primes. Theorem [2] implies 
that there is even a positive density of primes for which the conclusions of these 
three conjectures are false: 

Corollary 1 We have 



5 3>1 (1, 12) = ~ <f 3>1 (5, 12) = - 8 3>1 (7, 12) = - and 5 3)1 (11, 12) = 0. 
6 4 4 



Furthermore, we have 



Mi, io) 



i 



4,i (3, 10) 



1 



, 4,i (7, 10) 



1 



and 5 5) i (9, 10) 



1 



12 



4 



4 



12 



5 



In particular, the relative density of the primes for which the conclusion in Con- 
jectures 1.1-1. 4 fail are, respectively, 

MjMg) =n Mi, 12) = i Wg) = i MMo)J 

5(3) ' 5(3) 4' 5(5) 8' 5(5) 8' 

After Fermat various authors considered primes in arithmetic progressions divid- 
ing S a; b. Thus Sierpihski [T7] proved that every prime p = ±3(mod 8) divides 5*2,1 
and, furthermore, that no prime p = 7(mod 8) divides 5*2,1 • This result easily fol- 
lows on using that (-) = (— 1)^ _1 )/ 8 . Sierpihski states that M.A. Makowski has 
proved that infinitely many primes p = l(mod 8) divide 52,i (namely Makowski 
notices that the prime factors of the numbers of the form 2 2 " + 1 with n > 3 have 
the required property) and ends his paper with stating the problem of whether 
there are infinitely many primes p = l(mod 8) not dividing S^i- Subsequently, 
using results on the biquadratic and octavic residue character of 2, this problem 
has been independently resolved by A. Aigner pQ and A. Brauer [2]. Brauer 
shows for example that the infinitely many primes p = 9 (mod 16) which can be 
represented as 65x 2 + 256xy + 256y 2 all do not divide 5 2j i (the number of such 
primes < x is of order 0(x / \f\ogx) by a result of G. Pall [14], and thus this set 
has natural density zero). Using the first entry of Table 6 we infer that there 
many more primes not dividing 5*2,1: l/6th of all primes p = l(mod 8) do not 
divide S^i- 



3 The density written as infinite sum 

In order to evaluate 5 a> b{c, d) we will make use of the following result. 

Theorem 3 Let a, b, c, d be positive integers with c > 1 and d > 1 coprime. Let 
a c denote the automorphism of Q(Cd) determined by o- c ((d) = Q- The density 
o~a,b(c,d) of primes p = c(mod d) such thatp\S a ^ exists and satisfies 

MM) = g(p^-p^), (i) 

where 

r{j) = {\ lf =<• an d, similarly, r'(j) = { 1 lf a M = ld - 

y otherwise, ' 10 otherwise. 

Furthermore, Theorem^ holds true with 5 a> b{c,d) as given by (QJ). 

Proof. In case ord p (r) is defined we can define the index, i p (r), as (p— l)/ord p (r). 
Note that it equals [F* : (r)]. There is a unique j > 1 such that 2 J_1 | \i p {r). Let 
Pj denote the set of primes p such that 2 J_1 | \i p (r). Note that U^ 1 Pj equals, with 
finitely many exceptions, the set of all primes and that the are disjoint sets. 
Now note that for a prime p in Pj we have that ord p (r) is even if and only if 
p = l(mod 2 J ). Thus, except for finitely many primes, the set of prime divisors 
of 5 Qj fe satisfying p = c(mod d) is of the form W^LiQj, where 

Qj '■= {p '■ p = c(mod d),p = l(mod 2 J ), p e Pj}. 
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It is an easy observation that n\i p (r) if and only if p splits completely in Q(Cn, r l ^ n ). 
Using this observation and writing 's.c.' below to mean that the prime is split 
completely, we infer that 

Qj = {p : p = c(mod d),p s.c. in Q(C2i, 1 ),but not s.c. in Q(C2j r 1 / 23 )}. 

On invoking the Chebotarev density theorem, it is then found that the set Qj 
has a natural density that is given by 

HQ,)- T(3) Af> 



[N, : Q] [ATj : Q] ' 

On proceeding as in the proof of Lemma 8 of [15J it is then found that for 
ah < log 2//3 x and [d, 2 J ] < y : = log 1//6 xj log log x, and any number A > 0, we 
have 

Q j (x) = 5(Q j )U(x) + A (-^ r -). (2) 

Mog x' 

Thus 

3>X [d,2i]<y [d&]>y 

where n(x;m,n) denotes the number of primes p < x such that p = n(mod m) 
and Cj is any integer such that Cj = c(mod d) and Cj = l(mod 2 J ) if such an 
integer exists and 1 otherwise. A minor modification of the proof of Lemma 2 of 
[7] then yields that 



cd) W = £ ft(x) + (^£). ( 3 ) 

[d,2i]<» l0g X 



Using Lemma [2] we find that 

00 1 oA 

E *«M = o(2»£ |pp) = o( 7 ). (4) 

On combining (J2]), (j3J) and (jlj), the result is then obtained with 5 a t,(c, d) = 

ET=iWi)- D 

Remark 1. The algebraic side of the approach above (originating in Moree [7]) 
is not the traditional one to study the divisiblity of sequences S a ,b, but is chosen 
since it turns out to be easier to explicitly work out. The traditional approach 
rests on the observation that if p = 1 + 2- 7 (mod 2 J+1 ) for some j (which is uniquely 
determined), then ord p (r) is odd if and only if 

r (p-i)/a* = l(mod p), that is if 
and only if p splits completely in Q(C2j ', t 1 ^), see e.g. [12J for a sketch of the 
traditional approach. Note that (p — 1)/2 J ' is the largest odd divisor of p — 1 and 
so ordp(r) is odd if and only if ord p (r) divides (p — l)/2- ? . 

Remark 2. On GRH the existence of 5 a ,ft(c, d) was established by Moree [HI Theo- 
rem 1]. He showed under GRH that the set of primes p such that p = ai(mod d\) 
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and ordp(r) = a 2 (mod d 2 ) has a density S r (a\, d\\ a 2 \ d 2 ) and gave an expression 
for it in terms of field degrees and Galois intersection coefficients (r(j) and r'(j) 
in Theorem [3] are examples of such coefficients). Since S aj b(c,d) = S r (c, d; 0, 2), 
where r = a/b, it follows that 5 a ,b{c,d) exists under GRH. 

From our tables it is seen that 5 a ^{c, d) is always rational. Below a conceptual 
explanation for this is given. 

Proposition 2 The density 5 a ^{c, d) is always a rational number. 

Proof. We show that the sum in (JTJ) always yields a rational number. Note that 
Kj C Kj + i and K'- C K', +1 and hence the fields \im j^^Kj, lim j^^K^ exist. 
Denote these limits by K, K' . Note that K = K' . It follows that there exists jo 
such that t(J) = t'(J) and Kj = Kj = K = K' for every j > jo- By Lemma 
|2] it follows that there exist constants C\ and c 2 such that [Nj : Q] = Ci4 J and 
[TVj : Q] = for every j large enough. It follows that the terms with j large 
enough in ([T]) are in geometric progression and sum to a rational number. The 
terms are all rational and so 5 a ,&(c, d) is itself rational. □ 



4 Preliminaries on field degrees and field inter- 
sections 

The following facts from elementary algebraic number theory, for further details 
we refer to e.g. Moree [9], will be used freely in the sequel: 

1) a quadratic field K C Q(( n ) iff the discriminant of K divides n. 

2) Let Q(v / A) Q Q(Cn) be a quadratic fields of discriminant A and b be an integer 
with (b, n) = 1. Then Ofcl^^ = id. iff (j) = 1, with (7) the Krnecker symbol. 

In order to use Theorem [3] to compute 5 0) &(c, d), we first compute the degrees 
of the fields Nj, Nj for j > 1. This can be done directly or by using the general 
formula from Lemma 1 of [8] quoted below: 

Lemma 1 Put n t = [2 V2 ^ +1 , D(r )}. We have 



'fct,r v *) : Q] = ,, , where e(kt, k) - , , ... . , , 

7 J e(kt,k)(k,h) v ; \1 tfn t \kt. 

Using the lemma or otherwise, we compute the degrees of 

Nj = Q(C2„r 1 / 2J ",Cd) = QCC^C^r 1 /*" 1 ); 

to be as given in Lemma [2J The degrees turn out to be dependent on the follow- 
ing property which we call Cj : 

The property {Cj) holds if and only if D'\d', 5 < max(j, 5). 
Note that if D'\d', then (Cj) can fail only for finitely many j's. 
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Lemma 2 The degrees of Nj = Q((v,r 1/2j \( d ) and N- = Q(( 2 j,r 1/2 \( d ) over 
Q are given by: 



1 



<p(d) 



2 max(i,«5)-l ifj<\ + l; 

[Nj : Q] = { 2 max ^+i- x - 3 if j > \ + l and (Cj) holds; 

2 max(j,<5)+j-A-2 > a + 1, and (Cj) fails, 



( 2 max(j,«5)-l -y ■ < A . 

— -pr [N'j : Q] = < 2 max ^' <5 )+^ A - 2 if j > X and (Q) fto/ds; 
^ J ( 2 max ^^)+^ A - 1 if j > X and (Cj) fails. 

Remark 3. Equivalent form of (Cj). 

It will also be convenient to use the following version of (Cj) later. 

Property (Cj) holds if and only if either D(r )\d or D(r )\2 l d, D(ro) \ 2 l ~ l d for 
some I > 1 and j > I + 5. 

Equivalently, property (Cj) fails if, and only if, either D(tq) \ 2 l d VI > or 
D(r )\2 l d, D(r ) \ 2 l ~ x d for some I > 1 and j < I + 5. 

In the remainder of this section we assume that Q(y / r^) ^ Q(\/2). The case 
Q(v^o) 7^ Q(\/2) requires modification due to the ramification of 2 in cyclotomic 
extensions generated by large 2-power roots of unity and is discussed in Sections 
7 and 8. 

We need to determine precisely the set of all j > 1 for which r(j) = 1 and 
those for which r'(j) = 1. To this end we first determine the degrees of Kj,K'- 
over Q. 

Lemma 3 When 5 > 0, the degrees of Kj,K'j are given by the expressions : 

2 min(i,5) ifj<X + l; 
[Kj :Q} = { 2 min ^ 5 ) if j > X + 1 and (Cj) holds; 

2 min(i,5)-i l fj > \ + i a nd does not hold, 

2 min(j,<5) if ■ < X . 

[K'j :Q] = { 2 min ^) ifj>X and (Cj) holds; 

2 min(j,«5)-i l fj > \ an d (Cj) does not hold. 



Proof. When j < A + l, clearly r 1 ^ 23 1 is rational and, therefore, Kj = Q(C2 min W.*)] 
Similarly, K'- = Q(^ 2 ^u,s)) if j < ^- Further, note that Kj C K'- for all j. 
Writing Lj = Q(C 2 ,,r 1/2 ^ 1 ), and L'- = Q(( 2j , r 1 / 23 ), we have Nj = LjQ(( d ) and 
Kj = Lj n Q(( d ). Therefore, 

^ ;Q]= [^:Q][Q(0):Q] 



[Nj : Q] 

Similarly, N'j = L'^Q) and K'j = L) n Q(Cj). So, 



[*}=<»= [JVJ: 
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Using the above degree computations for Nj,N'j etc., we obtain the asserted 
expressions. 



For 5 = 0, the above formula has to be modified as we have used 4>(2 & ) = 2 s 1 . 
In this case, we get : 



Lemma 4 When 5 = 0, we have 
[Kj : Q] = 



2 if j > A + 1 and (Cj) holds; 
1 if either j < A + 1 or j > A + 1 and (Cj) fails, 



and 



2 if j > A and (Cj) holds; 

1 if either j < A or j > A and (Cj) fails. 



Remark 4. Since Kj is a subfield of Kj, it follows from the above degree compu- 
tation that Kj = in all cases except possibly when j = A + 1. For j = A + 1, 
we have Q(C 2 mm(A+M)) = K x+1 and the degree of K' x+1 over K x+1 is 2 if D'\d! 
and <5 < max(A + 1,5). If this latter condition (C\ +1 ) does not hold, then 
K\ +1 = K' x+1 . In other words, we have the following property : 

K j = K' j , r(j)=r'(j) V + 

We would like to actually write the fields Kj, Kj in a convenient form so that we 
can determine how the automorphism Q i— > ^ acts on them. Note that clearly 
the field Q(C 2 mm(j,i)) is always contained in Kj,K'- and its degree is either the 
whole or half of that of Kj,K'j . We look for a subfield of the form Q(( 2 ™™U,s)) 
or Q(C 2 minO',<s) , y/v) which has the full degree and will, therefore, have to be the 
whole field. 



Lemma 5 For j < \, Kj = Kj = Q(£ 2 mino\<s) ) . 
Furthermore, K\+i = Q(( 2 ™m^+l,s)) ■ 
For j >\ + \,Kj = K'j. 
For j > A + 1, K'j is : 

(a) Q(C 2 min(j,5)) if either D' \ d' or if 5 > max(j, 5) ; 

(b) Q(C2mino,o,v^o) ifD(r )\d ; 

(c) Q(C 2 mm(j,5), y/— r ) if D'\d' , 5 < 5o < max(j, 5), where r = u/v and 2 \ uv; 

(d) Q(C 2 minu,i), \Jl\i=i(^-)Pi) if D'\d',5 < 5 < max(j,5), where r = u/v with 
uv = 2 nLi Pi an d Pi > % f or i = 1, ■ ■ ■ , k. 

Proof. We know that Kj = Kj = Q(C2 min 0.*)) if either j < A or j > A + 1 and (Cj) 

fails. Also, K\ + i = Q(( 2 min(A + l,5) ) = K' x+1 unless (Cx+i) fails. In other words, we 

have to determine K'j only for those j > A for which (Cj) holds. 

Recall that the truth of (Cj) is equivalent to the property : 

either D(r )\d or D(r )\2 l d, D(r ) \ 2 l ~ x d for some 1 < / < 3 and j > I + 5. 

We examine each case separately. 

When D(r )\d, we have y/r^ G Q(Cz) and so, y/r\> G Kj. 

Moreover, if 5 > 1, then [Q(C 2 mi„o, 5) , ^/ro) : Q] = 2 mh < j '^ = [Kj : Q], except in 
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the case when Q(y/r~o) = <Q>(v2) which we have excluded in this section. Also, 
when 5 = 0, [Q(0=d) : Q] = 2 = \K' d : Q]. Therefore K' } = Q(C 2 »o-tf.*), V^o) if 
D(r )\d. 

When D(r )\2 l d, D(r ) \ 2 l ~ 1 d for some 1 < I < 3 and j > I + 8, it means that 
D'\d' , do — S + 1. If r = u/v, note that Q(^/ro) = Q(y/uv). Now, if uv is odd, 
it has to be = 3 (mod 4) since otherwise D(r ) = uv which cannot divide 2 l d 
without dividing d. Also then D(r Q ) = Auv = AD' , D'\d', 5q = 2 = 5 + I means 
that I = 1 = 5 or I = 2, 5 = 0. In case uv = 3(mod 4), we have a/— r € Q(v / ^) 
as the discriminant of Q(v /— r o) — — wi> = -D' which divides d' and hence divides 
d. Therefore K' d = Q(C 2 ^u,s), V^), when D(r )\2 l d, D(r ) \ 2 l ~ l d for some 
1 < I < 3 and j > I + 6 and r = u/v with -ut> odd. Here, we have used the fact 
that since j > 5o = 2, £ 4 (and hence y/—r~o) belongs to L'-. 

When uv = 2s with s > 1 odd, then D(r Q ) = Auv = 8s , 5 = 3, D' = s . Also 
5 = 8 — I = 3 — I and s = D'\d'. Thus, if s = Yl i==1 Pi, then \fi e Q(C Pl ... Pfc ) C 
Q(Crf); where i := rii=i(~)Pi- We have used the fact that \[2, i G Q(Cs) an d that 
j > So = 3. Hence K'j = Q(C 2 ^u,s), Vi) when uv is even and D(ro)\2 l d, D(r ) \ 
2 l ~ 1 d for some 1 < I < 3 and j > Z + 5. □ 

An immediate consequence of the previous lemma is the following result on the 
values of r(j) and r'(j). 

Lemma 6 If j < A + 1, Zaen r(j) = 1 min(j, 5) < 7. 

//j > A + 1 and «/ either D' \ d' or 5q > max(j, 5), then r(j) = 1 min(j, 5) < 7. 
J/ j > A + 1 and D(r )\d, then r(j) = 1 min(j, 5) < 7 and (^feel) = 1. 
If j > A + 1 and D'\d', 5 < 5o < j with uv odd where ro = u/v, then r(j) = 1 
min(j, 5) < 7 and ( D ^~ r °^ ) = I. 

If j > A + 1 and D'\d' , 5 < 5q < j wift w = 2 JI^Li where r = u/v and p^ 's 

odd primes, then r(j) = 1 <^ min(j, 5) < 7 and ( ^-1^ _ ^ 
VFe nai>e r'(j) = r(J) for j ^ A + 1. 

J/ e^ner £>' f cf or o" > max(A + 1,5), Znen r'(A + 1) = 1 min(A + 1,5) < 7. 
IfD(r )\d, then r'(A + 1) = 1 ^ min(A + 1,5) < 7 and (^) = 1. 
If D'\d',5 < 5 < A + 1 with uv odd where r = u/v, then r'(A + 1) = 1 4=> 
min(A + 1, 5) < 7 and (^=^) = 1. 

If D'\d', 5 < 5o < A + 1 wift uv = 2 \^ =l Pi where r = u/u , i/ien r'(A + 1) = 1 

• n _l 1 aw w . g(nii(j^) s . 
mm(A + 1, 0) < 7 ana ( ^— ) = 1. 

5 Tables for the density S a ^(c, d) when Q(y/ro) ^ 
Q(V2) 

Recall that the density 5 a> b(c, d) is given by (CQ). Since the primes considered are 
in (f)(d) residue classes, it is more natural to compute the sum 

S:=^)M^)=^)E(^%|-^). (5) 
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Note that S gives the relative density of divisibility of S a ,b, that is 

= #{p < x : p = c(mod rf), p|ffa, b } 
a:->oo ^£{p < x : p = c(mod d)} 

Putting in the degrees of Nj, iVj we can simplify the sum in (jSj) as follows. 

Since [Nj : Q] = [Nj : Q] and r(j) = r'(j) for j < A, the terms corresponding 
to j < A do not contribute. Also r(J) = r'(j) for j > A+l, but r(A + 1) and 
r'(A + 1) may be different (only) when (Cx+i) holds. Therefore, we have : 

= r(A + i)2 1-max ( A+1 ' <5 ) — r'(A + i)2 1 - max ( A+1 > <5 ) 
+ 2 A+1 E,>a + i,(C,) fails r(j)2—(^)^' if (C A+1 ) holds 

+ 2 A+2 E 1 >a + uc,) holds r(j)2— ^)^' 



s = r (A + i)2l- max (*+l> 5 ) - T (A + l)2- max ( A + 1 ' 5 ) 

+ 2 A+1 E^a+i,^) fails r{j)2-™«**>-i if (C A+1 ) fails 

+ 2 A+2 E J >A + i,(c,) holds rtf)2-™*tf.0-i 



As the degrees of the fields Nj,Nj and the values of r'(j)'s depend on the 
following three conditions, is convenient to have 3 tables depending on them. The 
three conditions are : 

(A) D> f d>; 

(B) D'\d',5 <5; 

(C) D'\d',5 > 5. 

Let us first work out the expression for S in case A. 



Case A: D' \ d' 

Here, every (Cj) fails. In particular, 



— = r (A + i)2- max ( A + 1 '' 5 ) + 2 A+1 r ( J )2- max 0'' 5 )-J. 



0(2 5 ) 

Moreover, since Kj = Kj = Q(( 2 ^inu,s)) f° r all j > A + 1, we have : 
For all j > A + 1, r(j) = t'(j) and this is 1 if and only if min(j, 5) < 7. 

Thus, S = <j>(2 5 )2^ Ej>A 1 min(j,.5)<7 2—^ = 0(2 5 )2 A + 1 (S 1 + S 2 ), 

where Si is the sum over j < 5 and S2 is the sum over j > S + 1. 
We get 

5 X = V 2-^andS 2 = f^>max(A + i,W) rJ if5<7; 

A + l<^in( 7 , 5 ) I ° therWiS 

From this, it is easy to obtain Table 1. 
Case B : D'\d', 5 < 5 
Note that (Cj) holds for all j. 

Here K x+ \ = Q(C 2 min(A+i,«)) and lf A+1 = Q(C 2 min(A+i,«), y/ro)- 
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For all j > A + 1, we have Kj = Kj = Q(£ 2 mm(j,<s), a/^o)- 
Therefore, r(A + 1) = 1 if and only if min(A + 1,5) < 7; 
r'(A + 1) = 1 if and only if min(A + 1, 5) < 7 and (^1) = 1. 
Moreover, for j > A + 1, we have 

r(j) = r'(j) which is 1 if and only if min(j, 5) < 7 and (^p^) = 1. 

Hence, we have 

s 

+ 2 A+2 E J >A + i^0')2' max( ^ ) - J , 
which can be written down more explicitly as S = <p(2 s )(t 1 + t 2 + S ), where 



J_ _ r ^ _)_ -g 2 i-max(A+i,<5) - r'(A + i)2 1 ~ max ( A+1 ' <5 ) 



ti = 



2 l-max(A+l,5) j f min ( A + ^ <J) < 7; 

otherwise, 
_ 2 i-max(A+i,tf) if min ( A + 1, 5) < 7 and (^M) = 1; 
otherwise, 

oA+2 n-max{j,S)-j ■£ ( D(rp) \ _ i. 

z Z^j>A+l,min(j,5)<7 Z 11 I c / ~~ L ' 

otherwise. 
Further, So = S 01 + S 02 , where S i is the subsum where j varies over j < 5 and 
S'02 is the subsum where j varies over j > 5. We find 



Soi — 
and that 



2 a+2-5( 2 -i-a _ 2 -min(7,«)) if (M) = 1 and A + 2 < min(5, 7); 
otherwise, 



•>02 



2 A+2-2max(A+l !( 5) /g jf (^o) ) = 1 a nd 5 < 7 ; 

otherwise. 
From this, we obtain Table 2. 

Finally, we work out the expression for S in case C. We write r = u/v and 
t = — r or rii=i(^")Pi according as to whether uv is odd or uv = 2 fl^LiPi- We 
also write D(t) for the discriminant of the quadratic field Q(\/F). 

Case C : D'\d', 5 > 5 

Notice that there are finitely many j's for which the property (Cj) may fail in 
this case. Now 

7^ 0(V \ TS> ( Q(( 2 min(X+l,5)) ifA + l<5 ; 

K, +1 = Q(C 2min( , + M)), K x+1 = | Q(C2mm(A+M)) ^ otherwise . 
For all j > A + 1, we have 

if = = J < Q(C2»»-o.*)) if 3 <5 ; 
j j \ Q(( 2 minu,s), Vi) otherwise. 

So, we have r(A + 1) = 1 if and only if min(A + 1,5) < 7 and furthermore we 
have 

^ fmin(A+ 1,5) < 7 , A + 1< 5 ; 
^ ^ \ min(A + 1, 5) < 7, A + 1 > 5 , and = 1. 
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Moreover, for j > A + 1 with j < 8q, we have r(j) = r'(j) which is 1 if and only if 
min(j, 5) < 7. On the other hand, for j > A + 1 with j > 8q, we have r(j) = r'(j) 
which is 1 if and only if min(j, 5) < 7 and (^i) = 1. 
Therefore, we get S = (fi(2 s )(ti + t 2 + Si + S2), where 

h = r(A + l)2 1 - max(A+1 ' 5) ; 

-r'(A + i)2 1 - max ( A + 1 ' 5 ) if A + 1 > S ; 
-r(A + i)2- max ( A+1 '< 5 ) if A + 1 < 5 ; 

S, = 2 A+1 J2{2-™^> s )-i : j > A + 1, j5 , min(j, 5) < 7}; 

S 2 = 2 x+2 J2{2~ m&x{j ' S) ~ j :j>X + l,j> 8 ,m.m(j,5) < 7} if (^) = 1 and, is 
0, otherwise. 

Putting in the values of r(A + 1) and r'(A + 1), we obtain 



2 l-max(A+l,5) if min ^ + ^ g) < T , 

otherwise, 



_ 2 i-max(A+i,5) a \ + 1 > 5 , min(A + 1,5) < 7, (^) = 1; 
t 2 = <j _ 2 -max(A+i,5) if A + 1 < 5 , min(A + 1, 5) < 7; 
otherwise. 
Finally, as before, we break up each of Si and £2 into two subsums over j < 5, 
respectively, over j > 6. So, we have Si = Su + 5*12, where 

S n = 2 X+1 - S ■ min (7, <*) > j > A + 1}; 

s = f 2 A+1 £{4~ J ' ■S > j> max(A + 2, 5 + 1)} if 0" < 7; 
to otherwise. 

Similarly, we have S 2 = £21 + £22, where 

s 21 = 0, S22 = ( 2A+2 : j ^ max ( A + 2 ' *<>)} if 5 ^ i and (^) = 

1 otherwise. 
On evaluating these expressions further we obtain Table 3. 



6 The intersection fields when Q(y/rS) = Q(V^) 

Next we consider the case where r = 2 or 1/2. Note that the discriminant of 
Q(v / 2) is 8 and that a/2 belongs to the cyclotomic field Q(Cs) (indeed y/2 = 
Cs + Cs" 1 )- Also note that Q(*> v^2) = Q(Cs) (we have Cs = (* + 1)/V^)- For j > 1 
we consider as before the degrees of the fields Nj, Nj. The earlier expressions in 
Lemma [2] are valid and, in fact, simplify to give: 

Lemma 7 TTte degrees of Nj = Q(Ca», r 1 / 2 '" 1 , a™d = Q(C 2 ., r 1/2i , over 
Q are gruen fry ; 

i f 2 max (-?> 5 )- 1 ifj<\+ 1; 

— — [iVj : Q] = < 2 max ^' <5 ) + '''- A - 3 »/ j > A + 1 and 3 < max(j, 5); 
^ > { 2 max (^)+^'- A - 2 ifj > A + 1 and 3 > max(j, 5), 
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2 max( J ,<5)-l % j ■ < A; 



— — [JVj : Q] = < 2 max ^ 5 )+-?- A - 2 if j > X and 3 < max(j, <5); 
^ ^ [ 2 max ^ 5 )+^- A - 1 z/j > A and 3 > max(?, 5) . 

The fields Kj = Q(C 2 , , r 1 ^' 1 ) n Q(0) and = Q(C^, r 1 / 2 ') n Q(0) are to be 
determined. This is where the computation gives different values from Lemma [3j 
However, the method of evaluation is the same and the degrees turn out to be : 

For j > X + 1, 

r 2 2 if j < 2, 5 > 3; 

[Kj : Q] = < 2 min ^< 5 )- 1 if either j > 3, 5 > 1 or j < 3, 1 < 5 < 2; 

1 1 if 5 = 0. 

For j > A, 

r 2 J ' if j < 2, 5 > 3; 

[Kj : Q] = < 2 min ^' 5 )- 1 if either j > 3, 5 > 1 or j < 3, 1 < 5 < 2; 

1 1 if 5 = 0. 



As we have evidently, l£j = Q(C 2 mm(j,<5)) for j < A + 1 and for every j, y^jnfati.*), 
is a subfield of Kj, we have the following result : 

Lemma 8 We have Kj = Q(( 2 minu.s)) for all j unless X = 0, j = 2, 5 > 3. 

In £/ie exceptional cases X = 0,j = 2,5 > 3, we have K 2 = Q(C 2 ^u,s) , v2) = 

Q(2,v / 2)=Q(Cs)- 

Further, we have Kj = Q(( 2 minu,s)) /or a// j unless X < j < 2,5 > 3. T/ie 
exceptional cases here are : either X = 0,j = 1,5 > 3 or A < 1, j — 2, 5 > 3. VFe 
/md t/ie following intersection fields: 

f A = 0, j = 1, 5 > 3, #J = Q(C 2 mina,*) , v^) = Qtv^); 
\ A < 1, j = 2, 5 > 3, K' 2 = Q{( 2 ^u, s) ,V2) = Q(i, y/2) = 



7 Tables for the density when Q(v^o) = Q(v^) 

Let S be defined as in ([5]). We divide its computation into four cases : 

(A) 5 < 3; 

(B) 5 > 3 and A > 2, 

(C) 5 > 3 and A = 1, and 

(D) 5 > 3 and A = 0. 



Case A : 5 < 3 



Then Kj = Kj = Q(( 2 minu,s)) for all j. Thus r(j) = r'(j) for all j and, this is 1 if 
and only if min(j, 5) < 7. It turns out that S = 0(2 <5 )(t 1 + t 2 + t 3 ), with 

2 -ma X (A+i,5) if ,\ < 1, min(A + 1,5) < 7; , f 1/8 ifA = 0,5< 7 ; 

1 t 2 



otherwise, ' ~ 1 otherwise 

t 



( 2 A+2-2max(A+l,2) if £ < 7; 

\ otherwise, 
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where t\,t 2 , t% correspond, respectively, to the terms in (jHJ) with j = A + 1, 
A + 2 < j < 3, j > max(3, A + 2) and j > max(3, 5 + 1). From this, we obtain 
Table 4. 

Case B : 5 > 3, A > 2 

Once again, Kj = Kj = Q(( 2 ^u,s)) for all j. Note that (Cj) always holds true. 
We obtain 

min(j,<5)< 7 



where 

2 1-S _ 2 A+2-5-min( 7 ,5) \ + 2 < min( 7 , 5)] 

otherwise, 



t 2 

i5\ + ,„fnS\ 



2A+2-2max(A+l,<5) ^3 jf £ < ry- 

otherwise, 



with cp(2 )ti, ip(2 )t 2 the subsum over j < S, respectively j > 5. 
Case C : 5 > 3, A = 1 

Here, we need to observe that when S\d, the Galois automorphism Q \— > Q of 
*d) fixes a/2 if and only if c = ±l(mod 8). We obtain 

S __r(A + l) r'(A + 1) , oA+2 ^ r(j) , qA+2 ^ r(j) 



+ 2 A+2 y YL. + 2 A+2 y 



cz)f2 <5N ) 2 s - 1 2 5 ' 1 2 max (J' ,5 )+J ^ 2 max (J' 6 )+-?' ' 

3<j<<5 j>max(2,5) 

which can be written as t 1 + t 2 + 1 3 + t 4 say, where further evaluation yields that 

h= T2 1 - 5 if2< 7 ; ) h= f-2 1 - 5 if3< 7 ; 

\ otherwise; ' \ otherwise; 
r 2 ^_ 2 3-«W) if3 < 7; iand . = |23-V3 if5< 7 ; 



, ciiiu t 4 . 

otherwise; [ otherwise. 

Table 5 is obtained from cases B and C. 

Case D : 5 > 3, A = 

As in the previous case, we need the fact that when 8\d, the Galois automorphism 
Cd >-> Q of Q(Cd) fixes y/2 if and only if c = ±l(mod 8). 
We find that S = 0(2*) (ti + t 2 + h + h), where 

2 1 - 5 if c= ±3(mod 8); = f 2~ 5 if 3 < 7 ; 
otherwise, " \ otherwise, 

^ o- 5 - 2 2 - 5 - min ^' 5 ) if 3 < min( 7 , 5); dt f 2 2 ~ 2S /3 if 5 < 7 ; 

otherwise, 1 1 otherwise, 



where ti,t 2 ,ts,t^ correspond, respectively, to the terms in ([5]) with j = 1, j = 2, 
3 < j < 5 and j > max(3, 5 + 1). This yields us Table 6. 
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8 Extremal densities 



We have < (p(d)5 ajb (c, d) < 1. In this section we are interested when 5 a ^(c, d) = 
and when 5 a> b{c, d) = l/ip(d). The following elementary result shows that if 
l(mod (d,2 A+1 )), then 5 a , b (c, d) =0. 

Lemma 9 Ifp\(a,b) and p\S a ,b, then p = l(mod 2 A+1 ). 

Proof. For a prime p put r(p) = (p — l)/(p — 1, /i). If p \ (a,b) and p\ab, then 
P t "Sa,&, so we ma y assume that p \ oh. Since r T ^ = {r^f^ = l(modp) by 
Fermat's little theorem, it follows that ord p (r)|r(p). If p is to divide S a ^ then 
r(p) must be even and so z/ 2 (p — 1) > A + 1. □ 

Theorem 4 a) Suppose that S a> b(c, d) = 0. This happens if and only if 

i) A > 7 and 5 > 7; 

or 

A = 7 - 1, £ > 7, £>(r )|d and (^pl) = 1. 
Moreover, if 5 a: b(c, d) = 0, then there are at most finitely primes p = c(mod d) 
dividing the sequence S a $. 

b) Suppose that 5 a ^(c,d) = l/ip(d). This happens if and only if 

i) A = 0, 5 = 0, D(r )\d and (^1) = -1; 
or 

ii) mm(7, 5) > A, Z}(r )|d and ( D ^°- ) ) = —1. 

Moreover, if £ a b (c, d) = l/ip(d), then there are at most finitely primes p = 
c(mod d) noi dividing the sequence S a> b. 

Proof. For a prime p put r(p) = (p — l)/(p — 1, /i). The first parts of both (a) and 
(b) follow on inspection of the Tables. Let us prove the second part of (a) now. 
If A > 7 and 5 > 7, we claim that r(p) is odd. Indeed, writing p = c + qd, and 
c — 1 = 2 7 co with Co odd, we have p — 1 = 2 7 co + 2 s qd'. Therefore, ^(p — 1) =7 
since 5 > 7. Now, (p — l,h) = (p — 1, 2 A /i') which has 2-adic valuation 7 since 
A > 7. Therefore r(p) is odd in the case (i) of (a) of the theorem. Since clearly 
ord p (r)|r(p), it then follows that p \ S a ^. Finally suppose we are in case ii. 
Suppose that p > 2 is a prime satisfying p = c(mod d) and such that p does not 
divide oh. Then, by the properties of the Kronecker symbol, 

where the first symbol is the Legendre symbol and fo denotes the reduction of ro 
modulo p. It follows that 

r 2(p-l,h) ^ l( m0( J p) ? 

and so ord p (r)|r(p)/2. We claim that r(p)/2 is odd. Now p — 1 = 2 7 co + 2 5 gd' 
which has 2-adic valuation 7 because <5 > 7. On the other hand, 2(p — 1, h) = 
2(p- 1, 2 A /i') = 2(p- 1, 2 7 " 1 /i / ) which has 2-adic valuation 1 + (7 - 1) = 7. Thus, 
r(p)/2 is odd and so p f S a ,b- 
b) The proof is similar; let us consider (i) first. 
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As 5 = A = 0, we have h is odd and r = r$. If p > 2 is a prime not dividing afe, 

then 

ro = / £?(r ) \ = / D(r ) \ = _^ 

p V p J V c / 

by assumption. Thus, Tq P = — l(mod'p), which implies that r^ p_1 ^ 2 = 

— l(modp) and therefore, that p\S aj b- Finally suppose we are in case ii. Writ- 
ing p = c + qd, and c — 1 = 2 7 c with c odd, we have p — 1 = 2 7 c + 2 <5 gc? / . 
Therefore, v 2 (p — 1) > min(5, 7). Now, v 2 (p — l,h) — v 2 (p — 1, 2 A /i') = A, since 
r> 2 (p — 1) > min(7,5) > A. Therefore, we have that is odd while r(p) is 

even; that is, 2 (p-lh) ls a P os itive integer. Once again, we have for each prime 
not dividing lab that 

( 5,. ( 5fe ! ) . ( ^ ) ._ L 

Thus, (r^ 1)/2 )(^W = -l(modp). But then r 2 ^- 1 ^) = (r^ -1 ^ 2 )^ 1 ^ = 

— l(mod p), which means that pjS'a^. □ 

Example. 1) By case ii of (a) we infer that £3,1(11, 12) = (cf. Conjecture 1.1 of 
Fermat). 

2) By case ii of (b) we infer that <^(8)<52,i(±3, 8) = 1 (easily proved using (2/p) = 
(— ~~ 1 ^ 8 ), cf. the paper by Sierpifiski [T7j . 

Perhaps a more illuminating phrasing of the above theorem is the following. 
Theorem 5 For a prime p put r(p) = (p — l)/(p — 1, h) . 

a) We have 6 at b(c,d) = if and only if r(p) is o<id or 2||r(p) and (^) = 1, /or 
all but finitely many primes p = c(mod d) . 

b) We have 5 at b(c,d) = l/ip(d) if and only if for all but finitely many primes 
p = c(mod d) we have that r(p) is even and (y) = — 1. 

Conclusion: if the density is extremal, then this can always be explained by ele- 
mentary arguments not using more than quadratic reciprocity and, furthermore, 
the associated set of exceptional primes is at most finite. 

Remark 5 (uniform distribution) . It is generally not true that the primes dividing 
S a ,b are uniformly distributed over the residue classes modulo d. However, there 
are some cases where we have uniform distribution. For example, if d is odd and 
D{tq) \ d, then the primes in any residue class mod d which divide S a ,b have the 
same density. 



9 Some numerical experiments 

For each entry in Tables 1-6 an example with parameters a and 6=1 was choosen 
and below we give the value of <5 a ,i(c, d) according to the tables on the one hand, 
and an approximation to this that consists of the first six decimals of the ratio 

#{p <Pm-P = c(mod d), p\S aA } 
#{p < Pm ■ P = c(mod d)} 
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where p m denotes the mth prime and m = 2097152000 « 2 • 10 9 . As a rule of 
thumb an approximation of 5 a ,i{c,d) obtained in this way by looking for prime 
divisors amongs the primes should have an accuracy of about n(p m ; d, c) -1 / 2 . We 
clearly observed in our experiments that for larger d the accuracy tends to be 
less (and the same holds for the run time). 



Test cases for Table 1 



Residue class 


a 


<f)(d)5 a>1 (c, d) 


Experimental value 


17 mod 56 


3 2 


5/6 


0.833200 • • • 


17 mod 56 


3 8 


1/3 


0.333317- •• 


1 mod 21 


5 


2/3 


0.666592 • • • 


7 mod 20 


3 4 








7 mod 20 


3 3 


1/2 


0.500015- •• 



Test cases for Table 2 



Residue class 


a 


<j)(d)5 atl (c, d) 


Experimental value 


9 mod 28 


7 2 


1/3 


0.333312- •• 


5 mod 12 


3 2 


1 


1 


1 mod 15 


5 


1/3 


0.333257- •• 


7 mod 15 


5 


1 


1 


1 mod 12 


3 


2/3 


0.666657- •• 


5 mod 12 


3 


1 


1 


11 mod 20 


5 4 








13 mod 24 


3 


1/2 


0.500006 • • • 


13 mod 56 


7 


1 


1 


7 mod 20 


5 2 









Test cases for Table 3 



Residue class 


a 


<fr(d)5 a>1 (c, d) 


Experimental value 


1 mod 12 


6 


11/12 


0.916693- •• 


5 mod 12 


6 


3/4 


0.749989 • • • 


1 mod 12 


6 2 


5/6 


0.833362- •• 


5 mod 12 


6 2 


1/2 


0.499996- • • 


7 mod 12 


6 


1/2 


0.500038- •• 


11 mod 28 


14 2 








7 mod 12 


6 4 








7 mod 30 


6 2 


5/12 


0.416679- • • 


11 mod 30 


6 2 


1/4 


0.250055- • • 


7 mod 30 


6 4 


1/12 


0.083321- •• 


11 mod 30 


6 4 


1/4 


0.250055- •• 


7 mod 15 


6 


17/24 


0.708336 • • • 


11 mod 15 


6 


5/8 


0.624999- • • 


7 mod 15 


6 4 


1/12 


0.083321- •• 


11 mod 15 


6 4 


1/4 


0.250055- • • 
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Test cases for Table 4 



Residue class 


a 


<j>{d)5 a Ac, d) 


Experimental value 


5 mod 14 


2 


17/24 


0.708327-- ■ 


5 mod 12 


2 


11/12 


0.916652 • • • 


7 mod 12 


2 


1/2 


0.499961 • • • 


7 mod 12 


2 2 








5 mod 6 


2 2 


5/12 


0.416673- •• 


5 mod 12 


2 2 


5/6 


0.833331 ••• 


5 mod 6 


2 8 


1/24 


0.041672 • • • 


5 mod 12 


2 4 


1/6 


0.166685- •• 


7 mod 12 


2 4 









Test cases for Table 5 



Residue class 


a 


(j)(d)S aA (c, d) 


Experimental value 


5 mod 24 


2 4 








17 mod 24 


2 4 


1/3 


0.333372 • • • 


17 mod 48 


2 4 


2/3 


0.666740 • • • 


17 mod 96 


2 4 


1/2 


0.500145- •• 


41 mod 48 


2 4 








17 mod 24 


2 2 


2/3 


0.666659 • • • 


7 mod 24 


2 2 








5 mod 24 


2 2 


1 


1 


17 mod 32 


2 2 


3/4 


0.750049 • • • 



Test cases for Table 6 



Residue class 


a 


(j)(d)5 aA (c, d) 


Experimental value 


9 mod 40 


2 


5/6 


0.833411- •• 


7 mod 8 


2 








5 mod 8 


2 


1 


1 


9 mod 16 


2 


3/4 


0.749983- •• 
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